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Abstract 

Massive multiple-input multiple-output (MIMO) networks, where the base stations (BSs) are equipped 
with large number of antennas and serve a number of users simultaneously, are very promising, but 
suffer from pilot contamination. Despite its importance, delayed channel state information (CSI) due to 
user mobility, being another degrading factor, lacks investigation in the literature. Hence, we consider 
an uplink model, where each BS applies zero-forcing decoder, accounting for both effects, but with the 
focal point on the relative users’ movement with regard to the BS antennas. In this setting, analytical 
closed-form expressions for the sum-rate with finite number of BS antennas, and the asymptotic limits 
with infinite number of BS antennas epitomize the main contributions. In particular, the probability 
density function of the signal-to-interference-plus-noise ratio and the ergodic sum-rate are derived for 
any finite number of antennas. Insights of the impact of the arising Doppler shift due to user mobility 
into the low signal-to-noise ratio regime as well as the outage probability are obtained. Moreover, 
asymptotic analysis performance results in terms of infinitely increasing number of antennas, power. 
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and both numbers of antennas and users (while their ratio is fixed) are provided. The numerical results 
demonstrate the performance loss in various Doppler shifts. An interesting observation is that massive 
MIMO is favorable even in time-varying channel conditions. 
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Index Terms 


Delayed channels, multiuser multiple-input multiple-output (MIMO), massive MIMO, zero-forcing 
(ZF). 


I. Introduction 

The rapidly increasing demand for wireless connectivity and throughput is one of the mo¬ 
tivations for the continuous evolution of cellular networks 0, 0. Very large multiple-input 
multiple-output (MIMO) has been identified as a new promising breakthrough technique aiming 
at achieving higher area throughput in wireless networks 0-0. Its origin is found in Q, 
and it has been given many alternative names such as massive MIMO, hyper-MIMO, and full- 
dimension MIMO systems. In the typical envisioned architecture, each base station (BS) with an 
array of hundreds or even thousands antennas, exploiting the key idea of multi-user MIMO (MU- 
MIMO), serves tens or hundreds of single-antenna users simultaneously in the same frequency 
band, respectively, under coherent processing. This difference in the number of BS antennas N 
and the number of users K per cell provides unprecedented spatial degrees of freedom that leads 
to high signal gains, allowing at the same time low-complexity linear signal processing techniques 
and avoiding inter-user interference due to the (near) orthogonality between the channels. 

On a similar note, zero-forcing (ZF) processing is regarded as a low-complexity alternative of 
maximum-likelihood multiuser detector and “dirty paper coding” Q, especially, when the BSs 
are equipped with massive antenna arrays. A lot of research has been conducted on single-cell 
systems with ZF receivers Q, but the main current interest has shifted to practical multi-cell 
scenarios, where pilot contamination degrades the system performance 0, [j^. 

Despite that the theory of massive MIMO has been now well established (see Q and references 
therein), an important question that has been overlooked is how the performance of massive 
MIMO topology is affected by the relative movement of users. This scenario is of high practical 
importance, in urban environments, where users move rapidly within a geographical area. The 
main challenge in time-varying environments is to perform robust channel estimation, when the 
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propagation channel changes over time. The dynamic channel behavior was modeled in terms of a 


stationary ergodic Gauss-Markov block fading channel model 110|-[ 121, where an autoregressive 
model was combined with the Jakes’ autocorrelation function that captures the time variation of 
the channel. Motivated by the above observation, this paper explores the robustness offered by 
massive MIMO against the practical setting of user mobility that results to delayed and degraded 
channel state information (CSI) at the BS, and thus, imperfect CSI. Such consideration is notably 
important because it can provide the quantification of the performance loss in various Doppler 
shifts. 

A limited effort for studying the time variation of the channel due to the relative movement of 
users has been conducted in [10|, where the authors provided deterministic equivalents (DEs|^ 
for the maximal ratio combining (MRC) receiver in the uplink and the matched filter (MF) 


in the downlink transmission. Fortunately, this analysis was extended in [11|, [12| by deriving 
DFs for the minimum mean-square error (MMSF) receiver (uplink) and regularized zero-forcing 
(downlink) and by making a comparison regarding their performance. In this paper, we consider 
a generalized uplink massive MIMO system. Based on the aforementioned literature, we propose 
a tractable model that encompasses ZF receivers and describes the impact of user mobility in 
a cellular system with BSs having conventional and very large number of antennas, or even 
large number of both antennas and users, which stands in contrast to the previous works. The 
following are the main contributions of this paper: 


Contrary to [I5|, we consider more practical settings where the channel is imperfectly 
estimated at the BS. The effects of pilot contamination and time variation of the channels are 
taken in to account. The extension is not straightforward because apart of the development 
of the model, the mathematical manipulations are hampered. Apart of this, the results are 
contributory and novel. 

We derive the probability density function (PDF) of the signal-to-interference-plus-noise 
ratio (SINK) and the corresponding ergodic sum-rate for any finite number of antennas 
in closed forms. For the sake of completeness, the link of these results with previous 


*The deterministic equivalents are deterministic tight approximations of functionals of random matrices of finite size. Note 
that these approximations are asymptotically accurate as the matrix dimensions grow to infinity, but can be precise for small 
dimensions. 
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known results is mentioned. Furthermore, a simpler and more traetable lower bound for the 
aehievable uplink rate is derived. 

• We elaborate on the low signal-to-noise ratio (SNR) regime, in order to get additional 
insights into the impact of Doppler shift. In particular, we study the behaviors of the 
minimum normalized energy per information bit to reliably convey any positive rate and 
the wideband slope. 

• We present a simple expression for the outage probability, being an important metric in 
quasi-static models. 

• We investigate the asymptotic performance presented by very large MIMO (N —)■ cx)) as 
well as large MIMO in terms of DEs (iV, K —)■ oo). This analysis aims at providing accurate 
approximation results that replace the need for lengthy Monte Carlo simulations. 

Note that, although all the results incur significant mathematical challenges, they can be easily 
evaluated. Nevertheless, the purpose of DEs is to provide the deterministic tight approximations, 
in order to avoid lengthy Monte-Carlo simulations. 

Notation: We use boldface lowercase and uppercase letters to denote vectors and matrices, 
respectively. The notation (■)^ stands for the conjugate transpose, and || ■ || denotes the Euclidean 
norm of a vector, while denotes the pseudo-inverse of a matrix. We use the notation x ^ y 
to imply that x and y have the same distribution. Einally, we use z ~ CM (0, E) to denote a 
circularly symmetric complex Gaussian vector z with zero mean and covariance matrix S. 


II. System Model 

Consider a cellular network which has L cells. Each cell includes one iV-antenna BS and K 
single-antenna users. We consider the uplink transmission. The model is based on the assumptions 
that: i) N > K, and ii) all users in L cells share the same time-frequency resource. Moreover, 
we assume that the channels are frequency flat and they vary from symbol to symbol, while 


during the symbol period they are considered constant due to the channel aging impact [10| (we 
will discuss about the channel aging model later). The channel vector giik[n] E between 

the kth user in the i\h cell and the (th BS at the nth symbol undergoes independent small-scale 
fading and large-scale fading. More precisely, giik[n] is modelled as 


OiM = \/Jhkhiik[n], 


( 1 ) 
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where I5iik represents large-seale fading, and huk G ~ CM (0, Iat) is the small-seale fading 
veetor between the Zth BS and the kth user in the zth eell. 

Let G be the zero-mean stoehastie data signal veetor of K users in the ith eell 

at time instanee n (pr > 0 is the average transmitted power of eaeh user, and Xi[n] ~ CM (0, Iat)). 
Then, the x 1 reeeived signal veetor at the Zth BS is 

L 

VM = y/M'^Gii[n]Xi[n] + Zi[n], l = l,2,...,L, (2) 

i=l 


where Gii[n] = [gmin],... ^ denotes the ehannel matrix between the K users 

in the ith eell and the Ith BS, and Zi[n\ ~ CA/^(0,lAr) is additive white Gaussian noise (AWGN) 
veetor at the Zth BS. 

To eoherently deteet the signals transmitted from the K users in the /th eell, the BS needs CSI 
knowledge. Conventionally, the /th BS ean estimate the ehannel via uplink training. We assume 


that the ehannel remains eonstant during the training phase [10|. In general, this assumption is 
not praetieal, but it yields a simple model whieh enables us to analyze the system performanee 
and to obtain initial insights on the impaet of ehannel aging. Furthermore, the impaet of ehannel 
aging ean be absorbed in the ehannel estimation error. 

During the training phase, in eaeh eell, K users are assigned K orthogonal pilot sequenees of 
length T symbols {r > K). Owing to the limitation of the eoherenee interval, the pilot sequenees 
have to be reused from eell to eell. We assume that all L eells use the same set of orthogonal 
pilot sequenees. As a result, the pilot eontamination oeeurs 0. i- Denote by ^ G 
(r > K), be the pilot matrix transmitted from the K users in eaeh eell, where the kih. row of 
’i' is the pilot sequenee assigned for the k\h user. The matrix ^ satisfies Then, the 

N X T reeeived pilot signal at BS I is 


rtrr 


n = 


+ Zf[n] 


/ — 1, 2,..., L, 


(3) 


i=l 


where the subseript O*'' implies the uplink training stage, ptr — ^Pr, and Zf[n] G is 

spatially AWGN at BS I during the training phase. We assume the the elements of Zf[n] are 
i.i.d. CM (0,1) random variables (RVs). The MMSE ehannel estimate of gukin] is given by 

1 


9iik[n\ =MkQik y^Sijkin] + ^ 

\i=l VPtr 



( 4 ) 
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where Qik^ + Aifcj and ~C7V(0, Iat) represents the additive noise. 

With MMSE ehannel estimation, the ehannel estimate and the ehannel estimation error are 
uneorrelated. Thus, giik[n] ean be rewritten as: 


duA^] =9iikH+9iik[ni 

where~ CA/" (O, [(3iik - Auk] liv) andp,^fc[n] ~ CAT (O, ), with= 


(5) 


^lik 


E, = l Ajfe + l/Ptr ’ 


are the independent ehannel estimate and ehannel estimation error, respeetively. Note that I3iik, 
Alik, and Qik are independent of n V/, i, and k. 

Besides pilot eontamination, in any eommon propagation seenario, a relative movement takes 
plaee between the antennas and the seatterers that degrades more ehannel’s performanee. Under 
these eireumstanees, the ehannel is time-varying and needs to be modeled by the famous Gauss- 
Markov bloek fading model, whieh is basieally an autoregressive model of eertain order that 
ineorporate two-dimensional isotropie seattering (Jakes model). More speeifieally, our analysis 
aehieves to relate the eurrent ehannel state with its past samples. For the sake of analytieal 


simplieity, we eonsider the following simplified autoregressive model of order 1 [10| 

9iik[n] = agiik[n - 1] + eiik[n], 


( 6 ) 


where gukin — 1] and eiik[n\ ~ CAf {0, (1 — Ajfcliv) are uneorrelated, denoting the ehannel at 
the previous symbol duration and the stationary Gaussian ehannel error veetor due to the time 
variation of the ehannel, respeetively. Note that a = Jq {27rfoTg), where Jo(-) is the zeroth-order 
Bessel funetion of the first kind, fo and are the maximum Doppler shift and the ehannel 
sampling period. Basieally, a, whieh is assumed known at the BS, eorresponds to the temporal 
eorrelation parameter that deseribes the isotropie seattering aeeording to the Jakes’ model. In 
partieular, the maximum Doppler shift equals /d = ^, where v (in m/s) is the relative 
veloeity of the user, c = 3 x 10®m/s is the speed of light, and fc is the earner frequeney. 

Substituting Q into we obtain a model whieh eombines both effeets of ehannel estimation 
error and ehannel aging as follows: 


9Hk[n] = agiikin - 1] + eiik[n] 

= agiik[n-l]+eiik[n\, (7) 

where gukin - 1] and eiik[n\ = agukin - 1] +eiik[n\ ~ CAf (^0, (^^Hk - liv) are mutually 

independent. More eoneretely, we define Gii[n\ = [giii[n], ■ ■ ■ ,giiK[n]] G and Eu = 
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7fc = - 


a^Pr 


a'^Pr 


E. 


L 




+ 



2 • 


( 12 ) 


as the combined channel matrices from all users in cell i to BS 1. 


[eiii[n],.. .,eiiK[n]] G C 

In particular, Gii[n] can be expressed as [15|: 

Gii[n] = Gii[n]Di, 

where ZJ. = diag{|j,|j,...,|f}. 

Making use of Q, we can rewrite the received signal yi[n] at the Zth BS (Z G [1, L\) as 

L L 

yi [n] = ay/p^'^Gii [n - l]Xi [n] + y/j^’^Eu [n]Xi [n] +Zi [n\. 


( 8 ) 


(9) 


Moreover, we assume that the Zth BS uses the ZF technique to detect the signals transmitted from 
K users in its cells. With ZF, the received signal vector yjn] is pre-multiplied with a~^G\i [n-1], 
where G\i[n — 1] is a Ff x iV matrix representing the pseudo-inverse of Gu[n — 1]: 

L 

riln] = y/p~rXi[n\ -f - l]Gii[n - l]Xi[n] 




-t-a ^^/j^y^Gll[n-l]Eli[n]Xi[n]+a ^G]i[n-l]zi[n], 


( 10 ) 


2=1 


Then, the fcth element of r; [n] is used to detect the signal transmitted from the kth user. The 
post-processed received signal corresponding to the kth user is 


rik [n] = y/VrXik N + \Gii [n-1] Gii[n-l]Xi[n] 


L 

2=1 


Eu[n\xi[n\ + - 

k OL . 


Gii[n-l] zi[n], 


( 11 ) 


where the notation [A\^ refers to the fcth row of matrix A, and is the kth element of Xi[n\, 

i.e., it is the transmit signal from the Zcth user in the /th cell at the nth time slot. Treating ( [TT] ) 
as a single-input single-output (SISO) system, we obtain the SINK of the transmission from the 
kth. user in the Zth cell to its BS in ( [T^ shown at the top of the page. The SINK is obtained 
under the assumption that the /th BS knows the denominator value of ( fT^ . This assumption is 
reasonable since this value is just a scalar (which can be estimated). 



























III. Achievable Uplink Sum Rate 


In this section, we provide the sum-rate analysis for finite and infinite number of BS antennas 
taking into account the aforementioned effects. 


A. Finite-N Analysis 

Proposition 1: The uplink SINK of transmission between the kth user in the /th cell to its 
BS, under the delayed channels, is distributed as 

d a^PrXk[n — 1] 


7fc 


(13) 


a^PrCXk[n - 1 ] + PrYkin] + 1 ’ 

where X^ and are independent random variables (RVs) whose PDFs are, respectively, given 
by 


PVfc (t) = ■ 




N-K 


, X > 0 


{N-K)\j3m VPukJ 
py, {y) =J2 2 /^ 0 ’ 


p=l q=l 


(14) 

(15) 


where Ak = diag Dn, G C 


••KLxKL 


D, 


kk 


with Dii & KxK diagonal matrix having elements 
= ( f^iik — oPj^iik), as well as g {Ak) denotes the numbers of distinct diagonal elements 


of Ak- Similarly, Pk,i,Fk, 2 , ■■■■, Fk,Q{Au) the associated distinct diagonal elements in decreasing 
order and (Alfc) are the multiplicities of pk,p, while Xp^q {Ak) is the (p, g)th characteristic 

coefficients of Ak, as defined in Definition 4]. Regarding (7, it is a deterministic constant: 

A (I3iik\‘^ 

[ 0uk ) ■ 

2 


Proof: Division of each term of ( [T^ by g\i [n- 1] 


leads to 


7fc = - 




. k 


-2 


a^PrC 


-2 


Gu[n-l] _ +PrEi=i 


YM 


(16) 


+ 1 


A 

\g^ \n 11 

G,-\n 11 

0 

iO“X and Y \n] ^ 


-‘-J 

k 




Since 


Gl[n-l] 


Gu[n - l]Gu[n - 1] 


)- 


si 

1 

_1 



kk 


k 


-2 


has an Erlang dis¬ 


tribution with shape parameter N — K + 1 and scale parameter j3iik [17|. Then, 

-2 


yt 


Gii[n-l] 


'Xk[n - 1]. 


(17) 
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r=0 


n+r 


(n + r)* 6”+'’ ^ (n + r)"'”*''’^ a6 ^ 


,s=0 


Q^n+r+lQ^m—n—r 


e-' -;^l-+^-l ^! Y^+r-s-l^^^n+r 

^ ^ a*a™-'*(a/a + l) 


—S—W—1 


a 


s+l 


. ( 20 ) 


Furthermore, eonditioned on 


Gii[n — 1] ,Yi[n] is a zero-mean eomplex Gaussian veetor with 


eovarianee matrix Du whieh is independent of 




independent of [n — 1] . Henee, X]i=i 

L i k 


G\i[n-l] Thus, Yi[n] ~ CM 


is the sum of KL statistieally independent 


Y,[n] 

but not neeessarily identieally distributed exponential RVs. Aeeording to Theorem 2], we 
obtain 

L 


^ ~ycN- 


(18) 


i=l 


Combining (fT^-(fT8l), we deduce ( fT3] ). ■ 

Remark 1: In the general case, the PDF of the uplink SINK ( [T3] ) accounts for both the effects 
of pilot contamination and Doppler shift. More specifically, the time variation of the channel 
decreases both the desired and interference signal powers by a factor of with comparison to 
the zeroth Doppler shift case, thus, degrading the SINK. 

Remark 2: Increasing the relative velocity, the SINK presents ripples with peak and zero 
points following the behaviour of the Jo(-) Bessel function. In the marginal case of a = 1, i.e., 
when there is no relative movement of the user, ( fT3] ) expresses the downgrade of the system 
only to pilot contamination. Especially, if we assume very long training intervals (orthogonal 
pilot sequences) and no time variation, which is not practical in common scenarios with large 


number of antennas and moving users, our result coincides with [13, Eq. (6)]. At the other end, 
high velocity meaning a —)■ 0 leads to zero SINK. 

Corollary 1: Consider the high uplink power regime. We have 

MXk[n - 1 ] 


d 


MCXk[n - 1] + TfcN’ 


as Pr —)■ oo. 


(19) 


This corollary brings an important insight on the system performance, when increases 
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N-K 




-ePuk<-'^Pr(C + -i) 


t=0 

N-K-t 

^ E , 

( Pllka'^Pr {C + 1) 

N-K r 


hlkOl^ {C + 1) f^llkO^'^Pr {C + 1) Pk,p PiikOt^ {C + 1)^ 

1 


(u 1)!( Pll_^_^ r{g)U[g,g+l + N-K-t-u, 


f^k^pPr 


( 22 ) 






t=0 
N-K-t 


1 


1 


fiiikoi^C (3iika‘^PrC Pk,p jSiika'^C. 


- ^^'^^N-K-t-u T{q)U (q,q + l + N -K-t-u, ^ 

“=1 ( hlkOp'PrC 


P'k^pPr 


(23) 


asymptotically. As seen in ( [T9| ), there is a finite SINK eeiling due to the simultaneous increment 
of the desired signal power as well as of the interferenee and channel estimation error powers. 


Having obtained the PDF of the SINK, and by defining the funetion J'rn,n (a, b, a) as in ( |20| ) 
shown at the top of the previous page, where Ei (■) denotes the exponential integral funetion 


119 Eq. (8.211.1)], we first obtain the exaet Rik {pr,C() and a simpler lower bound Rl (Pr,®) 
as follows: 

Theorem 1: The uplink ergodie aehievable rate of transmission between the fcth user in the 
/th eell to its BS for any finite number of antennas, under delayed channels, is given by 


SiAkRpiAk) -y ( A \ l,-,rr o 

D / \ ^ ^ 

Rik{Pr, a) =2^ “^ 2 ), 


( 21 ) 


where Xi and X 2 are given by ( |22| ) and ( [23] ) shown at the top of next page, and where (/(•,•,•) 
is the eonfluent hypergeometrie funetion of the seeond kind p9| Eq. (9.210.2)]. 

Proof: See Appendix ■ 

In the case that all diagonal elements of Ak are distinct, we have q {Ak) = KL, Tp {Ak) = 1, 
and Ap^i {Ak) = Y\^=i q^p • The uplink rate beeomes 


KL N-K 


n 


KL 


k‘k 


Rit(p.,a) E E (Ar_/ 4 '_t)!(_i) 


’Srf'Pe^e 


N-K-t 


Pk,p 


(T 1 -X 2 ), 


(24) 


where Xi and X 2 are given by ( [25] ) and ( |26l ) shown at the top of next page. Note that, we have 
used the identity U (1,&, c) = e^x^~^r {b — l,x) [29 Eq. (07.33.03.0014.01)] to obtain ([24[). 
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N-K 

t=0 
N-K-t 


ehk<-‘^Pr(C+l) 


E 


e“-i)!(-i) 


f^llkOi'^ (C + 1) i^llkOt'^Pr {C + 1) Pk,p {C + 1) 

1 


=1 {(5uka^C + l) 


\N+l-K-t-n, 


l^k^pPr 


(25) 


N-K 




-e'Puk^'^prC 


-K-t 


t=0 '- 
N-K-t 

E , 

l^iikOt^C 


j3iikOi'^C (3iika^PrC f^k,p jSiikO.'^C. 


(«-l)!(-l) 


[N + l-K-t-u, 


P'k,pPr 


■ (26) 


-Pout (O^th) 


1, if 7th > 1/C 


ei^k)'rp{Ak)N-K t 


)^E' E E EOXp,Mk)§ir{s + qpuk{a^-a^C^,^ 

P =1 q^l s =0 ^ 


s+g 


if 7th < l/C. 


(29) 


Proposition 2: The uplink ergodic rate from the kth user in the /th eell to its BS, considering 
delayed channels, can be presented by a certain lower bound RL{pr,a): 


Rlk {Pr,Oi)>RL {Pr,a) 

( 


= logs 




1 + 


67 + 


{N-K)a^f}uk 


(^E E(/3z*fc-a%fc) + ^ 


(27) 


Proof: See Appendix ■ 

1) Outage Probability: Bearing in mind that we investigate a block fading model, the study 
of the outage probability is of crucial interest. Basically, it defines the probability that the 
instantaneous SINK 7 ^. falls below a given threshold value 7 th: 


Pout (Tth) (7fc E Tth) • 


(28) 
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Theorem 2: The outage probability of transmission from the kih. user in the Zth eell to its BS 
is given by ( |29l ), shown at the top of the page. 

Proof: See Appendix ■ 

B. Characterization in the Low-SNR Regime 

Even though Theorem renders possible the exaet derivation of the uplink sum-rate, it appears 
defieient to provide an insightful dependenee on the various parameters sueh as the number of BS 
antennas and the transmit power. On that aeeount, the study of the low power eomerstone, i.e., the 
low-SNR regime, is of great signifieanee. There is no reason to eonsider the high-SNR regime, 
beeause in this regime an important metrie sueh as the high-SNR slope S^o = linip^_ 5 .o 
[ (2^ is zero due to the finite sum-rate, as shown in ( [T9l ). 

1) Low-SNR Regime: In ease of low-SNR, it is possible to represent the rate by means of 
seeond-order Taylor approximation as 


. .. 

Rik (Pr, a) = Rik (0, a) Pr + Rik (0, a) ^ + o (pi) , 


(30) 


where Rik {pr ,«) and Rik {pr, a) denote the first and seeond derivatives of Rik {pr ,«) with 
respeet to SNR pr. In faet, these parameters enable us to examine the energy effieieney in the 
regime of low-SNR by means of two key element parameters, namely the minimum transmit 
energy per information bit, 

Pr 


and the wideband slope Sq |24|. Espeeially, we have 

1 


= lim 


.^Omin Pr—^0 Rik {^pry Cl) Rik (0, o) 


5n = 


Rik (0, a) 


^ln2. 


(31) 


(32) 


Rik (0, a) 

Theorem 3: In the low-SNR regime, the uplink sum-rate betweem the fcth user in the /th eell 
to its BS in a multi-eell system, assuming delayed ehannels, ean be eaptured by the minimum 
transmit energy per information bit, , and the wideband slope Sq, respeetively, expressed 

'Omin 

by 

Eh ln2 


A^Omin aHN-K+l)^iik 

-2{N-K+1) / {N-K+2) 


5n = 


a-‘+2a-C(lV-ii'+3) + ^ E T. 

p=l q=l !■ 


(33) 

(34) 


Proof: See Appendix 
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C. Large Antenna Limit Analysis 


In this section, we consider the large system limit by accounting for specific assumptions. 
Assuming constant transmit power and : i) the number of BS antennas N grows infinitely 
large, while K is fixed, ii) both the number of uses K and BS antennas N increase asymptotically 
by keeping their ratio k = ^ fixed; and with scaling the power with the number of antennas 
N, we obtain the corresponding SINRs, in order to scrutinize their properties. The purpose of 
this analysis is to exploit the reduction of the interference and thermal noise due to the property 
of orthogonal channels vectors between the BS and the users as iV —)■ cx) as well as to achieve 
increase of the sum-rate due to its dependence of N. 

i j iV — )■ oo with fixed Pr and K: Keeping in mind that an Erlang distribution with shape 
and scale parameters given by iV — A' + 1 and fiuk, respectively, can be related with the sum of 
independent normal RVs Wfin — 1], W2[n — 1],..., W2(N-K+i)[n — 1], Xfin — 1] as follows: 


W|n - 1] = ^ E W^\n - 1], 


2{N-K + 1) 


(35) 


i=l 


By the substitution of ( [35] ) into ( [T3| ) as well as by the use of the law of large numbers, the 
nominator and the first term of the denominator in ( [T3] ) converge almost surely to a^prf3iik/2 
and a‘^prCfiiik/2 as N ^ oo, while the second term of the denominator goes to 0. As a result, 
the deterministic equivalent of the SINR, 7 ^, when N ^ 00 , is expressed as: 


a.s. _ 

7fc 7fc 


1 


as —)■ 00 . 


(36) 


The bounded SINR is expected because it is already known that as the number of BS antennas 
tends to infinity, both the intra-cell interference and noise are cancelled out, while the inter-cell 
interference due to pilot contamination remains. 

2) K, N ^ 00 with fixed pr and n = N/K: In practice, the number of serving users K in 
each cell of next generation systems is not much less than the number of BS antennas N. In 
such case, the application of the law of numbers does not stand because the channel vectors 
between the BS and the users are not anymore parrwisely orthogonal. This in turn induces new 
properties at the scenario under study, which are going to be revealed after the following analysis. 
Basically, we are going to derive the deterministic approximation 7 ^ of the SINR 7 ^ such that 



7fc - Ik 


N —^00 


(37) 
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where 


a.s. 


N^oo 


> denotes almost sure eonvergenee. 


Theorem 4: The deterministie equivalent 7 ^ of the uplink SINK between user k and its BS 
with ZF deeoder is given by 


Ik = 


oi^hik (k - 1 ) 


a'^CPiik (k - 1) + Ylf=i :^Tr£) 


(38) 


Proof: Knowing that Yfn] ~ CM we ean write it as: 

where ai~ CM if), \k)- By substituting ( |35] ) as well as ( [39l ) into ( [13] ), we have 


Ik = 


Mpr^ Mi[n - 1 ] 


2 (iv-r!'+i) 


(39) 


(40) 


MprC^ Y, W^[n-l]+prY o-fDaUi + 1 


2 = 1 2=1 

Next, if we divide both the nominator and denominator of (|4^ by 2 (iV — iT + 1) and by using 


110 Lemma 1], under the assumption that Du has uniformly bounded speetral norm with respeet 


to K, we arrive at the desired result (38). 


Remark 3: Interestingly, in eontrast to ( |3^ , the SINK is now affeeted by intra-eell interferenee 
as well as inter-eell interferenee and it is independent of the transmit power. In faet, the former 
justifies the latter, sinee both the desired and interferenee signals are ehanged by the same faetor, 
if eaeh user ehanges its power. Note that the interferenee terms remain beeause they depend on 
both N and K; however, the dependenee of thermal noise only from N makes it vanish. As 


expeeted, ( [38] ) eoineides with ( |42| ), if S> iT, i.e., when k —>■ 00 , the SINK goes asymptotieally 
to XjC. 

Next, the deterministie equivalent sum-rate ean be obtained by means of the dominated 


eonvergenee [25| and the eontinuous mapping theorem |26| as 


Rik{Pr,a) -log2(l + 7fc) 


a.s. 


N^oo 


> 0 . 


(41) 


3) Power-Scaling Law: Let eonsider p^ = E/y/N, where E is fixed regardless of N. Given 
that Mk depends on ptr = we have from (|3^ that for fixed K and N ^ 00 , 


Ik 




Ilk 


(42) 


aME^C(3l, + V 

whieh is a non-zero eonstant. This implies that, we ean reduee the transmit power proportionally 
to I/a/ZV, while remaining a given quality-of-serviee. In the ease where the BS has perfeet CSI 
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Fig. 1. Sum spectral efficiency versus SNR for different N{a — 0.1 and a = 0.9). 


and where there is no relative movement of the users, the result (42) is identical with the result 


in [13|. 


IV. Numerical Results 

In this section, we present numerical results to verify our analysis by considering a cellular 
network with L = 7 cells and K = IQ users per cell. The coherence interval is T = 200 symbols 
(which corresponds to a coherence bandwidth of 200 kHz and a coherence time of 1 ms) and 
the length of training duration is t = K symbols. Regarding the large-scale coefficients I3iik, we 
assume a simple scenario: I3iik = 1 and fiuk = a,ior k = 1,... ,K, and i ^l. Note that a can be 
considered as an intercell interference factor. In all examples, we choose a = 0.1. Furthermore, 
we define SNR = pr. 

In the following, we will examine the sum spectral efficiency which is defined as: 

K 

k=l 


( 43 ) 
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a = 0.1, SNR = 0 dB 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Temporal Correlation Parameter (a) 


Fig. 2. Sum spectral efficiency versus a for different N{a = 0.1 and SNR = OdB). 


where Rik {pr, «) is given by ( [^ . 

Figure shows the sum speetral effieieney versus the SNR for N = 20, 50, and 100, at the 
intereell interferenee faetor a = 0.1, and the temporal eorrelation parameter ct = 0.9. The “Exaet, 
Analysis” eurves are eomputed by using ( |2T] ), the “Exaet, Simulation” eurves are generated via 
( fT^ using Monte-Carlo simulations, while the “Bound” eurves are obtained by using the bound 
result given in Proposition The exaet mateh between the analytieal and simulated results 
validates our analysis. It ean be seen from the figure that, the proposed bound is very tight, 
espeeially for large antenna arrays. Furthermore, we ean see that, at high SNR, the sum speetral 
effieieney saturates. This is due to the faet that when SNR inereases, both the desired signal 
power and intereell interferenee power are inereased. To improve the system performanee, we 
ean use more antennas at the BS. At SNR = 5dB, the sum speetral effieieneies ean be inereased 
by the faetors of 2.5 or 5.5 if we inerease N from 20 to 50 or from 20 to 100, respeetively. 

Next, we examine the effeet of the temporal eorrelation parameter on the system performanee 
as well as the tightness of our proposed bound given in Proposition Figure presents the 
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Fig. 3. Transmit power required to achieve 1 bit/s/Hz per user versus N{a = 0.1, a = 0.7 and a = 0.9). 

sum spectral efficiency as a function of the temporal correlation parameter a, at SNR = OdB, 
for N = 20, 50, and 100. We can see that the system performance degrades significantly when 
the temporal correlation parameter decreases (or the time variation in the channel increases). 
Half of the spectral efficiency is reduced when a reduces from 1 to 0.6. Furthermore, at low a, 
using more antennas at the BS does not help improve the system performance much. Regarding 
the tightness of the proposed bound, we can see that the bound is very tight across the entire 
temporal correlation range. 

Figure shows the transmit power, p^, that is needed to reach 1 bit/s/Hz per user. Here, 
we choose a = 0.1 and a = 0.7 or 0.9. As expected, the required transmit power decreases 
significantly when we increase the number of BS antennas. By doubling the number of BS 
antennas, we can cut back the transmit power by approximately 1.5dB. This observation is in 
line with the results of 

To further verify our analysis on large antenna limits, we consider Figure]^ Figure]^ shows the 
sum spectral efficiency versus the number of base station antennas for different values of a, and 
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Fig. 4. Sum spectral efficiency versus N for different a. 


for two cases: the transmit power, pr, is fixed regardless of N, and the transmit power is sealed 
as Pr = I/VN. The “Limits” eurves are computed via the results obtained in Seetion III-C. 
As expeeted, as the number of the base station antennas increases, the sum speetral effieiencies 
eonverge to their limits. When the transmit power is fixed, the asymptotie performanee (as 
N —)■ oo) does not depend on the temporal eorrelation parameter. By eontrast, when the transmit 
power is scaled as the asymptotie performance depends on a. 

Finally, we eonsider the outage performanee versus SNR at = 100, for different temporal 
correlation parameters (a = 1, and 0.9), and for different threshold values ( 7 th = 2, and 3). 
See Figure We can see that, the outage probability strongly depends on a. At SNR = 0 dB, 
by redueing a form 1 to 0.9, the outage probability inereases from 7 x 10“® to 5 x 10“^, and 
from 3 X 10“^ to 5 x 10“^ for 7 th = 2, and 3, respectively. In addition, the outage probability 
significantly improves when the threshold values are slightly redueed. This is due to the fact 
that, with large antenna arrays, the ehannel hardening oeeurs, and henee, the SINK eoneentrates 
around its mean. As a results, by slightly redueing the threshold values, we ean obtain a very 
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Fig. 5. Outage probability versus SNR for different a and 7 tii (A'^ = 100). 


low outage probability. 


V. Conclusions 

In this paper, we characterized the uplink performance of a cellular network taking into account 
both the well-known pilot contamination and the unavoidable, but less studied, time variation. The 
latter effect, inherent in the vast majority of propagation scenarios, stems from the user mobility. 
Summarizing the main contributions of this work, new analytical closed-form expressions for the 
PDF of the SINK and the corresponding achievable sum-rate, that hold for any finite number of 
antennas, were derived. Moreover, a complete investigation of the low-SNR regime took place. 
Neveretheless, asymptotic expressions in the large antenna/user limit were also obtained, as well 
as the power-scaling law was studied. As a final point, numerical illustrations depicted how the 
time variation affects the performance in various Doppler shifts for finite and infinite number 
of antennas. Notably, the outcome is that large number of antennas should be preferred even in 
time-varying conditions. 
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Appendix 


A. Proof of Theorem 

The uplink ergodic rate from the /cth user in the Zth eell to its BS (in bits/s/Hz) is given by 

Pra‘^Xk[n - 1] 


Rik(pr: a) =Exfc,y^<jlog 2 (l-l 

’log2 (u 



Pra^CXk [n-1] +PrYk [n] +1 

Px, (a:) PYk (y) dxdy 


PrO^X 


>0 JO 

Q{Ak) Tp(Ak) 

h (g-l)!(iV-A^)!A 


PrO^Cx + PrV + 1 

R^P,q (A) Pkl log2 e 

‘TV-X+l 


X 



Inl 


0 ao 


Profx 


PrO^Cx+Pry+l 


— 'X — y 

dxdy 


Q{Ak) Tp(Ak) -y ( A \ ..“a o 
I'Afcj p iog 2 e 




X 




0 Jo 


Pr-y+l 


— X —y 

x^-^ehk yi~^e^k,p dxdy 


=Xi 

°° foo / ry‘^Cr\ ^ 

In 11 + — -- I x^~^e ^uk j/i-ig f'fe.p dxdy 



0 ao 


P 


- ; 

■y+i) 


— 

=X2 


Q{Ak) Tp(Ak) -y ( A \ ..“a 1„„ £) 

Xp^q [Ak) Pk,p ^^§2 6 ,x -r ^ 

: 77777; 77777777777 “ ^2j • 

p=l q= 


^ (g-l)!(iV-i7)!/3™ 


(44) 


(45) 


We first derive Xi by evaluating the integral over x. By using p9| Eq. (4.337.5)], we obtain 

N-K „oo 

11 = 57 / [-f{vX’‘-‘e-f^y>Ei{f{y)) 

/—n J 0 


i=0 
N-K-t 


J2 


U=1 


where f{y) 


PrV+l 


yilkPrA(C+l) 


. Using [23 


Lemma 1] and [21 


y^l ^e^k,pdy^ (46) 

Eq. (39)], we ean easily obtain Xi 


as given in ( [22] ). Similarly, we obtain X 2 as given in ( [23] ). Substitution of Xi and X 2 into ( |45] ) 
coneludes the proof. 


















B. Proof of Proposition 

By using Jensen’s inequality, we have 
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Rik (Pr, a) = E {log 2 (1 + 7 fc)} = E log 2 1 + 


1/7a 


> log2 1 + 


= Rl (Pr, Oi) ■ 


E{l/7fc}. 

To compute Rl (Pr, oi), we need to compute E { 1 / 7 ^.}. From ( [T^ , we have 


Ik 




CT —Ve 


2 = 1 


a^Pr 


E 


Gii[n-l] ^Eu[n] 


Gii[n-l] 


= C' + —E 


= C- 



f3*|n-l] 


n 



k 

i 


L K 


, i=l k=l 


[l3lik — 0'‘^^lik) + 


Pr 


L K 




Pr 


(47) 


(48) 


{N-K)a^Puk Vtrt-i 

In the third equality of ( [48] ), we have considered the independence between the two variables, 
while in the last equality, we have used the following result: 

2 ' 


E 


vt 


Gu[n-l] 


= 


Xk[n - 1] 


roc ^-x/0uk f X \ _ 1 

/o {N-K)\^lXj^k) "" “ {N-K)^ak 


(49) 


Note that we have used p9[ Eq. (3.326.2)] to obtain ( |49l ). Thus, the desired result ( [T7| ) is obtained 
from ( |47] ) and ( |48] ). 


C. Proof of Theorem 

Clearly, from ( fT3] ), 7 ^. < 1/C. Thus, if 7 th > 1/C, then Pout ( 7 th) = 1- Next, we consider the 
case where 7 th < 1/C. Taking the probability of the instantaneous SINR 7 fc, given by ( [13] ), we 
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Rlk {Pr 


a 


ln2 


a'^Xk[n - 1] - 1] + 2<jfc (1 + Pr<;kf (l + a^PrXk[n - 1] + Pr<;k)^ 

(aVr (C' + 1) X‘^[n - 1] +pYk[n\ + 1)^ {a'^p'^CX^ln - 1] + prYk[n\ + 1)"^ 


, (56) 


can determine the outage probability as 

a^PrXk 


Pn„t. = Pr 


Xk Y Pr^k “t“ 1 
Tth iPr^k Y 1) 


Pr Xfc < 


'0 


N-K 

= 1 — gPrTth ^ ^ 
t=0 ' 


(y^Pr - ItYLoYPrC 


< 7th 

Yk ) PYk{y)^y 


-y 

eT-th 
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n-k( ^ 

E ■ 


4=0 

e(Ak) Tp{Ak) N-K 


\7th 

TT 


= l_ep.4 ^ ^ ^Xp^q{Ak 
p=l t=o 


yY — ]pYkiy)dy 

Pr/ 

-q ( 

Yk,p \7tii J 

{q - 1) t\ 


X I y'^ f 2/ + ~ ) ^y 


Jo 

eiAk)rp{Ak)N-K t 

-7th ^ sr^sr^l t 


Pr 


t 


=i-e-E E EE(>..wyys. 

P=1 g=l f=0. = oW 


(50) 


where 7th = Puk — Q!^C'7th), and where in the third equality, we have used that the eumulative 
density funetion of (Erlang variable) is 


Fxi^ix) = Pr (Xfc < x) 


= 1 


exp 


X 


f^llk 


Y,(£) 


4=0 


(51) 


The last equality of ( [SO] ) was derived after applying the binomial expansion of {y Y I/prY 
and ^ Eq. (3.351.1)]. 


D. Proof of Theorem 

The initial step for the derivation of the minimum transmit energy per information bit is to 
eover the need for exaet expressions regarding the derivatives of Rik {pr,a). In partieular, this 
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can be given by 

Rik (pr, “) = 

ja'^Xk [n - l]/(aVrC'^fc [n - 1] +PrYk [n] +1)| 

(a*p=(C+l).Y2|n-l|+p,n|n] + l) j 

Easily, its value at = 0 is 

Rik (0, a) = — Ejsf^ {a^Xfc[n — 1]} • 


( 52 ) 


(53) 


Taking into account that [n — 1] is Erlang distributed, its expectation can be written as 


Ex, {Xk[n-l]} = {N-K + l)^Puk. 


(54) 


Substituting ( |54l ) and ( |5^ into ( [3T] ), we obtain the desired result. 

The second derivative of Rik {pr, a), needed for the evaluation of the wideband slope, is given 
by ( [5^ shown at the top of the previous page, where ca: = a‘^CXk[n — 1] + Yk[n\. Hence, 
Rik (0, a) can be expressed by 


Rik ( 0 , a) = {a^Xl[n - 1 ] + 2a^CXl[n - 1 ] 


+2a^Xk[n - l]Yk[n]] . 


(57) 


The moments of Xk [n — 1] are obtained by means of the corresponding derivatives of its moment 
generating function (MGE) at zero (0), i.e.. Ex, {XJ![n — 1]} = (0). Thus, having in 

mind that the MGE of the Erlang distribution is 


Mx, (*) 


1 


— Azfcf j 


N-K+l ’ 


we can obtain the required moments of Xk [n — 1] as 


(58) 


Ex, {X,2[n-l]} = Mg)(0) 

r(iV-iT + 3) ^, 
T{N-K + 1)^^^^ 

Ex, {X|[n-l]} = Mg(0) 

r(iV-iT + 4) ^3 
T {N - K + 


(59) 


( 60 ) 
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In addition, since Xk[n — 1] and Yk[n] are uneorrelated, we have {^k[n ~ = 

Exfc {Xk[n — 1]} In other words, it is neeessary to find the expeetation of Yk[n]. As 

aforementioned, the PDF of Yk[n] obeys ( [T5| ) and has expeetation given by definition as 


^YA^k[n]}= I ypYAy)^y 

JQ 

oXk) TpiAk) 




p=l q=l 

e(Ak) Tp(Ak) 

p=l q=l 


(?- 1)! Jo 
^kU 


y<]Ql^k,p dy 


(61) 


where we have used [19 Eq. (3.326.2)] as well as the identity r(g + l) = q\. As a result, 

,101. Finally, substitution of the 


Rik (0, a) follows by means of ( |59| ), 
yields the wideband slope. 


and (57) into (32) 
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